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Executive  Summary 


The  goal  of  this  DARPA  TIFT  project  was  to  demonstrate  an  efficient  room  temperature  source 
for  THz  imaging  systems  based  on  cascaded  optical  down-conversion  in  an  orientation-patterned 
GaAs  microstructure,  incorporated  into  an  optical  resonant  cavity,  and  pumped  by  a  compact 
solid-state  laser.  The  main  achievements  of  this  Project  are:  (1)  quasi-phase-matched  (QPM) 
GaAs  microstructures  with  periodic  reversal  of  crystalline  orientation,  including  orientation- 
patterned  GaAs  (OP-GaAs),  diffusion-bonded  GaAs  (DB-GaAs)  and  optically-contacted  GaAs 
(OC-GaAs)  suitable  for  THz  generation  were  designed  and  fabricated  with  the  length  exceeding 
2  cm,  thickness  of  more  than  1  mm  and  optical  absorption  coefficient  less  than  0.01  cm'1,  (2) 
efficient  resonantly  enhanced  THz  wave  generation  with  GaAs  inside  the  cavity  of  a  near¬ 
degenerate  doubly-resonant  optical  parametric  oscillator  was  demonstrated  in  the  whole  range 
0.7  -3.5  THz,  with  an  average  THz  power  of  1  mW  at  2.8  THz,  and  (3)  cascaded  optical-to-THz 
frequency  conversion  was  demonstrated,  which  opens  up  the  possibility  of  exceeding  the 
Manley-Rowe  conversion  efficiency  limit  that  is  getting  more  than  one  THz  photon  per  optical 
photon. 
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1.  Introduction 

There  is  a  colossal  potential  for  exploiting  terahertz  waves  (1  THz  =10  Hz,  k=300  pm)  in  the 
fields  extending  far  beyond  the  realms  of  their  traditional  use,  like  astronomy,  study  of  Earth’s 
atmosphere,  high-resolution  spectroscopy  and  plasma  diagnostics.  These  new  applications 
emerged  in  the  last  two  decades  and  encompass  ‘friendly  X-rays’  real-time  imaging  (THz 
radiation  experiences,  in  many  occasions,  much  smaller  scattering  than  the  optical,  and  thus  can 
penetrate  many  materials;  yet  the  photon  energy  is  too  small  to  do  any  harm  to  living  organisms), 
sensing  and  spectroscopic  imaging  by  means  of  rotational-vibrational  spectroscopy,  because  of 
extreme  richness  of  absorption  ‘fingerprints’  in  the  THz  range,  biomedical  imaging  (identifying 
cancers,  particularly  skin  cancer),  pharmaceutical  industry  (classifying  molecular  polymorphs), 
as  well  as  broadband  wireless  communication  [1],  Another  emerging  field  is  nonlinear 
interactions  of  THz  waves  with  matter  and  nonlinear  terahertz  spectroscopy  [2,3,4], 

Regardless  of  a  great  variety  of  techniques  for  generating  THz  radiation  [5],  none  of  them 
provides  a  THz  source  which  is  simultaneously  (i)  compact,  (ii)  highly  efficient,  (iii)  broadly 
tunable  (octave  or  more),  and  (iv)  works  at  room  temperature.  THz  sources  based  on  nonlinear 
optical  effects,  in  particular  those  using  y(2)  nonlinearity,  such  as  difference  frequency  generation 
(DFG),  optical  rectification,  and  optical  parametric  oscillation  (OPO)  are  very  appealing  because 
they  can  capitalize  on  swift  progress  of  diode-pumped  lasers,  in  particular  fiber  lasers,  which 
would  allow  making  a  THz  system  very  compact.  Another  advantage  is  power  scalability:  unlike 
photomixers  and  dipole  photoantennas  [5],  nonlinear  optical  sources  do  not  suffer  from 
saturation  of  THz  output  power  and  typically  show  quadratic  dependence  vs.  pump  power.  It 
should  be  mentioned  that  THz  radiation  can  also  be  produced  from  optical  radiation  using  y(3) 
nonlinearity:  intense  THz  pulses  were  generated  via  four-wave  rectification  process  in  laser- 
created  plasma  in  ambient  air  [6,7]. 
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In  this  paper,  we  briefly  review  different  techniques  of  THz  generation  using  ya>  processes, 
including  quasi-phase-matched  interactions,  make  comparison  of  linear  and  nonlinear  properties 
of  crystals  suitable  for  THz-wave  generation,  consider  theoretical  limits  of  optical-to-THz 
conversion,  and  review  optical  rectification  experiments  performed  with  periodically-inverted 
GaAs.  Finally,  we  will  describe  and  review  first  experimental  results  on  a  new  approach  to 
optical  THz-wave  generation.  This  approach  is  based  on  frequency  mixing  in  periodically 
inverted  GaAs,  between  two  resonantly-enhanced  optical  fields  inside  a  cavity.  We  will  also 
discuss  recycling  of  optical  photons  via  cascaded  process,  which  makes  it  possible  to 
substantially  surpass  Manley-Rowe  conversion  limit. 


2.  Optical  THz-wave  generation  in  nonlinear  crystals 

2.1.  Frequency  conversion  using  DFG  and  OPO 

Generation  of  THz  waves  by  difference  frequency  mixing  using  laser  radiation  and  a  nonlinear 
crystal  was  first  performed  by  Zemike  and  Berman  in  quartz  crystal  inside  the  cavity  of  a  ruby 
laser  as  early  as  in  1965  [8],  Later,  tunable  THz  generation  by  DFG  was  demonstrated  in  other 
nonlinear  optical  materials:  both  collinear  (forward  and  backward)  and  noncollinear  phase 
matching  was  used  to  generate  broadly  tunable,  20-190cm'1  (0.6-5. 7  THz),  output  in  LiNbCf  [9], 
THz  DFG  was  achieved  with  collinear  phase-matching  in  birefringent  crystals  of  ZGP  [10,1 1] 
and  GaSe  [12],  and  with  noncollinear  phase  matching  in  cubic  zinc  blende  structure  crystals  of 
GaAs  [13,14]  and  GaP  [15],  as  well  as  with  collinear  phase-matching  in  GaP  [16]. 

Guided-wave  interactions  between  optical  and  THz  waves  can  confine  radiation  to 
transverse  dimensions  on  the  order  of  the  THz  wavelength.  Since  DFG  THz  efficiency  depends 
on  pump  power  density,  it  can  be  enhanced  in  waveguides,  thus  lifting  the  tradeoff  between 
confinement  and  diffraction.  In  addition,  DFG  phase  matching  in  a  THz  waveguide  made  of 
nonlinear  material  can  be  achieved  with  optically  isotropic  nonlinear  materials  like  GaAs,  due  to 
the  waveguide  dispersion  added  to  the  bulk  dispersion  of  a  crystal.  Planar  GaAs  waveguides 
(170  and  75 -pm  thick)  were  used  in  ref.  [17]  to  generate  0.5  and  1.2  THz  radiation  with  ~  mW 
peak  power,  by  mixing  outputs  of  two  Q-switched  CO2  lasers. 

Alternatively,  broadly  tunable  (0.3-3  THz)  quasi-  monochromatic  THz  radiation  can  be 
produced  through  THz-wave  parametric  oscillation  in  bulk  LiNbCf  crystal  with  a  single  fixed- 
frequency  optical  pump  [18,19,20].  In  this  case,  phase  matching  is  achieved  through  noncollinear 
three-wave  interaction;  the  pump  and  the  OPO  resonating  wave  lie  in  the  near  infrared,  while  the 
outcoupled  non-resonant  wave  is  the  THz  beam. 


2.2.  Optical  rectification  (OR) 

In  optical  rectification,  THz  output  is  produced  in  a  nonlinear  medium  via  difference-frequency 
mixing  between  Fourier  components  of  the  same  optical  pulse.  First  demonstrated  with 
picosecond  pulses  in  LiNb03  [21,22],  and  ZnTe,  ZnSe,  CdS  and  quartz  [21]  crystals,  this 
technique  was  later  extended  to  femtosecond  laser  pulses  [23]  which  made  it  possible  to  generate 
much  broader  bandwidths  in  terahertz  region  [24],  Typically,  to  generate  broadband  THz 
transients  via  optical  rectification,  thin  (~lmm  or  less)  electro-optic  crystals  are  used,  because  of 
constraints,  associated  with  the  mismatch  between  optical  group  and  terahertz  phase  velocities, 
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as  well  as  high  absorption  in  conventional  crystals  (LiNb03,  ZnTe)  at  THz  frequencies.  Optical 
rectification  was  also  demonstrated  in  zinc -blende  GaAs  [25,26]  and  GaP  [27]  crystals,  which 
have  intrinsically  much  smaller  THz  absorption.  In  the  latter  case,  GaP  was  collinearly  pumped 
by  an  ultrafast  (x~210  fs)  Yb-doped  fiber  laser  operating  at  1.055  pm.  A  120-MHz-repetition- 
rate  pulse  train  of  single-cycle  THz  radiation  with  6.5  pW  average  power  was  generated  using  10 
W  of  laser  pump  [27].  Optical-to-THz  conversion  efficiencies  in  optical  rectification  are 
typically  in  the  range  10"6  -10"9  [5],  However,  recently  OR  conversion  efficiency  of  ~  5x1 0~4  was 
demonstrated  in  LiNbCF  by  tilting  the  intensity  front  of  a  femtosecond  (fs)  800  nm  pump  pulse 
to  match  the  phonon-polariton  phase  velocity  (at  THz  frequencies)  to  the  group  velocity  of  the 
pump  pulses  in  LiNbCF  crystal.  The  authors  generated  near  single-cycle  pulses  with  an  energy 
per  pulse  of  240  nJ  [28]  and  10  pJ  (100  pW  average  power)  [29],  In  the  latter  case,  energetic  (20 
mJ)  fs  laser  pulses  at  10  Hz  repetition  rate  were  used  as  pump. 


2.3.  Quasi-phase-matched  interactions  applied  to  THz  generation 

Crystals  with  periodic  variation  of  the  x(2)  nonlinear  coefficient,  can  solve  the  problem  of  k- 
vector  conservation  in  optical-to-THz  conversion  process,  via  quasi-phase  matching  (QPM).  In 
the  case  of  optical  rectification,  on  the  other  hand,  QPM  crystals  increase  effective  length  of 
collinear  interaction  by  eliminating  destructive  interference  due  to  velocity  mismatch  between 
optical  and  THz  waves.  Optical  rectification  in  a  QPM  crystal  was  first  demonstrated  in 
periodically-poled  lithium  niobate  (PPLN)  [30].  The  authors  attained  narrow-band  THz 
generation  by  using  femtosecond  pulses  at  800  nm  and  PPLN  crystal  with  different  QPM  periods 
and  achieved  ~  10'5  conversion  efficiency.  PPLN  was  cryogenically  cooled  (T=18  K)  to  reduce 
THz  absorption.  In  contrast  to  velocity-matched  single-cycle  OR,  optical  rectification  with  QPM 
crystals  gives  rise  to  a  multi-cycle  narrow-band  terahertz  radiation  with  the  center  frequency 

Qo=2nc/  AAn  (1) 

and  spectral  width  corresponding  to  the  QPM  acceptance  bandwidth 

AQthz=2tt:c/ZA«  (2) 

where  A  is  QPM  period  and  An=nmz-ng  opt  is  a  mismatch  between  the  THz  phase  (nm:)  and 
optical  group  ( ng  „pt)  indices.  Each  inverted  domain  of  a  QPM  nonlinear  crystal  contributes  a 
half-cycle  of  the  THz  pulse  and  thus  the  THz  wave  packet  has  as  many  oscillation  cycles  as  the 
number  of  QPM  periods  over  the  length  of  the  crystal  (Fig.l).  Also,  PPLN  was  used  in  surface 
emitting  geometries  for  THz  generation  using  both  OR  (fs  pulses)  [31]  and  DFG  (ps  pulses)  [32], 
Lately,  THz-wave  generation  was  demonstrated  in  periodically-inverted  GaP  [33]  and  GaAs 
[34,35],  In  fact,  III-V  semiconductors  are  very  attractive  for  QPM  THz-wave  generation  because 
of  several  appealing  properties,  namely  (i)  small  THz  absorption  coefficient  (smaller  by  an  order 
of  magnitude  than  in  commonly  used  EO  crystals:  lithium  niobate,  ZnTe,  CdTe),  (ii)  large 
coherence  length  due  to  small  mismatch  between  the  optical  group  and  THz  phase  velocities,  and 
(iii)  high  thermal  conductivity.  In  Ref.  33,  THz- waves  were  obtained  by  DFG  in  a  periodically- 
inverted  GaP  stack,  produced  by  direct-wafer-bonding  technique,  and  the  GaP  crystal  was 
pumped  with  10-nanosecond  pulses  near  1.55  pm. 


Thus,  optical  THz  generation  in  nonlinear  crystals  provides  a  variety  of  methods  to  produce 
THz  output  with  an  average  power  from  nW  to  mW,  and  peak  power  up  to  megawatts.  For 
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example,  single-cycle  THz  pulses  with  peak  power  of  5  MW  were  demonstrated  by  OR  in 
LiNbOs  [29],  and  sub-ns  pulses  with  2  MW  peak  power  (at  1  THz)  -  by  noncollinear  DFG  in 
GaAs,  using  250-ps  CO2  laser  pulses  [36]. 


3.  Theoretical  limits  of  optical-to-THz  conversion  in  QPM 
crystals 


As  in  any  x(2)  parametric  down-conversion  process,  optical-to-THz  conversion  efficiency  scales 
as  Q2,  and  as  Q3,  in  the  case  of  tightly  focused  beams  [37,38],  In  particular,  for  the  plane-wave 
DFG  with  long  pulses  [37]: 


dfHz  — 


2Q2d2ffL 2 

*0C  noptnTHz 


(3) 


where  t|Thz  is  optical-to-THz  conversion  efficiency  by  intensity,  Q  is  the  angular  terahertz 
frequency,  de ff  is  the  effective  nonlinear  optical  coefficient,  L  is  the  length  of  the  crystal,  Iopt  is 
the  optical  intensity,  80  is  the  permittivity  of  free  space,  c  is  the  speed  of  light,  nopt  is  the  optical 
refractive  index,  and  nT Hz  is  the  THz  refractive  index.  Because  of  Q2  (or  Q3)  term,  THz 
conversion  efficiency  is  intrinsically  low,  as  compared  to  the  case  of  generating  near-  or  mid- 
infrared  radiation  with  the  same  technique;  to  achieve  close  to  unity  photon  conversion 
efficiency  in  DFG  process  at  1-2  THz,  one  needs  optical  pump  intensity  on  the  order  of  1 
GW/cm2  for  both  GaAs  and  LiNbOs,  provided  that  the  interaction  length  is  close  to  the  inverse 
THz  absorption  coefficient.  Moreover,  very  large  Stokes  shift  imposes  limitations  on  maximum 
conversion  efficiency  which  can  be  achieved  even  if  100%  photon  conversion  is  reached. 


3. 1.  Long  pump  pulses 

At  a  given  pump  power,  the  maximum  THz  efficiency  in  a  collinear  DFG  can  be  achieved 
with  tightly  focused  beams  (approximately  confocal  with  respect  to  the  THz  wave).  In  this  case, 
for  sufficiently  long  pulses  (ns  to  CW),  optical-to-THz  conversion  efficiency  scales  [38]  as  t|Thz 
~  RoptG  (Jeff  /«0pt  )L  ~  Ropt  G  (Jeff  /n0 pt  )/otTHz,  were  Popt  is  the  optical  pump  power;  we  also 
assumed  that  the  optimal  L  is  on  the  order  of  1/a  thz,  where  a  thz  is  THz  power  absorption 
coefficient  (THz  absorption  usually  dominates  over  optical  absorption).  Thus,  we  can  introduce  a 
figure  of  merit  for  long  pulses  as: 

FOMi=  (Jeff2/n0pt2)/ctTHz  (4) 

3.2.  Quasi  phase-matched  OR  and  DFG  with  short  pulses 

Parametric  frequency  conversion  is  intensity-dependent  and  it  might  seem  that,  from  the 
viewpoint  of  achieving  high  optical-to-THz  efficiency,  OR  is  more  attractive  than  DFG,  because 
of  access  to  higher  peak  powers  with  fs  pulses.  It  appears,  that  in  a  QPM  case,  for  sufficiently 
short  (xG<l)  fs  optical  pulses,  THz  conversion  efficiency  in  OR  scales  with  the  pulse  fluence, 
not  intensity  [38],  The  physical  reason  is  that  in  a  QPM  crystal,  fs  pulse  interacts  coherently  with 
the  THz  wave  only  over  the  temporal  walk-off  length,  associated  with  the  difference  in  optical 
group  and  terahertz  phase  velocities,  /w~  ex/ An  «L  (here  c  is  the  speed  of  light,  x  -  pulse  length); 
hence,  shortening  the  pulse  increases  the  peak  intensity  but  simultaneously  decreases  /w. 
Furthermore,  in  the  case  of  mixing  two  ps  pulses  (xG»l),  with  pulse  durations  such  that  the 
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temporal  walk-off  length  /w  is  smaller  than  the  length  of  the  crystal,  again,  conversion  efficiency 
is  fluence-,  not  intensity-  dependent,  and  has  the  same  scaling  factor  as  in  the  case  of  OR  with 
femtosecond  pulses.  In  this  event,  for  both  OR  and  DFG  (with  optimally  focused  beams),  t|Thz 
scales  with  pulse  energy  Uopt  as  t|thz  ~  £/0ptO  (de ff  In opt  )/  An,  and  a  figure  of  merit  for  energy 
conversion  efficiency  for  short  pulses  appears  as: 

FOM2=  ( <ieff2/«opt2  y  An  (5) 

As  calculations  show  [38],  for  QPM  GaAs,  t|Thz  ~0.1%  per  pj  at  1-2  THz,  for  both  fs  and  ps 
pulses  with  optimized  pulsewidths. 

3.3.  Kerr-limited  interactions 

In  the  previous  example,  conversion  efficiency  did  not  depend  on  crystal  length,  because  we 
assumed  that  for  shorter  crystals  we  apply  stronger  focusing.  This  is  true  when  pump  intensities 
are  moderate  and  we  neglect  high-order  nonlinear  effects,  such  as  two-  and  three-photon 
absorption,  and  nonlinear  refraction.  Assuming  that  the  pump  wavelength  is  such  that  there  is  no 
two-photon  absorption,  and  three-  and  higher  photon-number  absorption  negligible,  we  conclude 
that  the  limiting  factor  is  the  nonlinear  refraction  (Kerr  effect).  To  avoid  self-focusing  and  self¬ 
phase  modulation,  pump  intensity  Iopt  should  be  kept  low  enough  to  have  nonlinear  phase  in  the 
optical  beam  center  to  be  Acp  =  27i/kopt-  n2  IopfL  <  1  rad.,  where  kop,  is  the  optical  wavelength,  and 

«2  is  the  index  of  nonlinear  refraction  (by  intensity).  In  this  case,  it  is  advantageous  to  use  longer 

2  2 

crystals  (T~l/a  thz);  maximum  THz  conversion  efficiency  scales  as  Aopt(r/eff  /«op t  «thz)T/«2  ~ 

2  2 

^oPt(^eff  /«0pt  «thz)/<2thz«2-  Then,  a  figure  of  merit  for  Kerr-limited  conversion  is: 

FOM3=  Aoptr/eff2 / (»oPt2«THz  CC THz  »2)  (6) . 


4.  Comparison  of  nonlinear  crystals  for  THz  generation 

4.1.  Linear  properties 

Many  nonlinear  optical  materials  are  relatively  transparent  at  THz  frequencies  below  lowest 
phonon  resonance  (e.g.  at  8  THz  for  GaAs  and  5  THz  for  LiNbCf)  and  their  absorption 
monotonically  decreases  with  decreasing  frequency.  Absorption  can  be  affected  by  free  carriers, 
but  in  high  resistivity  materials  (e.g.  GaAs),  the  carrier  density  can  be  made  extremely  low  and 
not  contribute  to  absorption.  Some  weak  features  at  low  frequencies  are  due  to  multiphonon 
processes  involving  longitudinal  and  transverse-acoustic  modes  or  optical  and  longitudinal 
phonons.  Fig.  2  shows  THz  power  absorption  coefficient  for  seven  nonlinear  crystals  at  room 
temperature.  Most  of  the  data  are  based  on  recent  measurements  [2,15,14,39,40,41],  One  can  see 
that  GaAs,  GaP  and  ZGP  (ZnGeP2)  show  best  THz  transparency  below  4  THz. 

4.2.  Nonlinear  properties 

For  THz  frequencies  well  below  phonon  resonances,  the  effective  nonlinear  coefficient  dy\  for 
DFG  and  OR  is  nearly  constant  and  is  derived  from  the  linear  electro-optic  coefficient,  r,  using 
the  relation  [42]:  d^\=  -  n4/(4r j^),  where  n  is  the  optical  refractive  index.  For  a  number  of 
semiconductor  materials,  dy\  is  lower  at  THz  frequencies  than  in  the  optical  range,  while  for 
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ferroelectric  oxides  like  LiNbCb  -  there  is  an  opposite  trend.  Fig.  3  shows  frequency  dependence 
of  the  nonlinear  coefficient  du  for  GaAs,  on  a  broader  scale  of  wavelengths  [17]. 


4.3.  Figures  of  merit 

Linear  and  nonlinear  optical  properties,  as  well  as  figures  of  merit  for  the  crystals,  which  are 
transparent  in  the  0-4  THz  range,  are  listed  in  Table  1  [14,  43-48],  In  some  crystals  (ZnTe,  GaP) 
index  mismatch,  An,  can  be  close  to  zero;  in  this  case,  to  calculate  FOM2,  we  took  the  walk-off 
length  /w~  ex/ An  to  be  equal  (for  x~100  fs)  to  the  physical  length  L~  1/aiHz ,  since  /w  can  not 
exceed  L.  LiNb03  crystal  appears  to  have  the  largest  nonlinear  optical  coefficient  in  the  THz 
range,  however  GaAs  has  much  smaller  THz  absorption;  as  a  result  their  figures  of  merit  FOMi 
are  comparable. 

5.  QPM  GaAs  samples 

Three  different  types  of  GaAs  samples  with  periodically-  inverted  crystalline  orientation  were 
used  for  QPM  THz  generation:  (i)  diffusion-bonded  GaAs  (DB-GaAs)  [49],  produced  by 
stacking  and  bonding  together  alternately  rotated  [110]  GaAs  plates;  wafer  fusion  in  this  case 
creates  a  monolithic  body  with  periodic  change  in  the  nonlinear  coefficient,  (ii)  orientation- 
patterned  GaAs  (OP-GaAs)  [50]  grown  by  a  combination  of  molecular  beam  epitaxy, 
photolithography,  and  hydride  vapor  phase  epitaxy  (HVPE),  where  periodic  inversions  of  the 
crystallographic  orientation  are  grown  into  the  material,  and  (iii)  optically-contacted  GaAs  (OC- 
GaAs)  [51],  which  was  similar  to  DB-GaAs  since  the  technology  involves  separate  plates  of 
GaAs  which  are  brought  together,  in  a  clean  room  environment,  with  180°  rotation  around  [110] 
direction  between  neighboring  plates;  unlike  DB-GaAs,  the  plates  were  not  heated  to  create  a 
uniform  crystalline  structure.  While  DB-GaAs  and  OC-GaAs  provide  large  (-1  cm2)  apertures, 
OP-GaAs  (Fig.  4)  has  the  advantage  of  more  reproducible  technology  and  allows  precise 
lithographic  definition  of  QPM  gratings;  however  it  suffers  from  limited  aperture  (thickness)  of 
0.4-1  mm.  Transmission  loss  at  optical  (k~2  pm)  frequencies  was  measured  to  be  ~  6%  per  cm 
in  DB-GaAs.  With  the  newly-grown  OC-GaAs  and  OP-GaAs  samples,  transmission  loss  was 
measured  to  be  <1  %  per  cm.  Parameters  of  GaAs  samples  used  in  experiments  are  listed  in 
Table  2. 


6.  THz  generation  in  GaAs  by  optical  rectification 

6.1.  OR  with  fs  pulses  from  a  Tm-fiber  laser  at  2  jum 

In  order  to  eliminate  two-photon-absorption  in  GaAs,  one  needs  to  use  optical  pumping  with 
wavelengths  X>  1.75  pm,  corresponding  to  the  photon  energy  which  is  less  than  half  of  the 
bandgap  of  GaAs.  In  the  OR  experiment,  OP-GaAs  crystals  were  pumped  by  a  compact 
femtosecond  fiber  laser  operating  near  2  pm  (Fig.  5).  An  all-fiber  laser  produced  120-fs  pulses  at 
100  MHz  repetition  rate  with  the  average  power  of  3  W  (30  nJ  pulse  energy)  at  the  wavelength 
of  1980  nm.  In  this  fiber  laser  design,  the  output  of  a  mode-locked  Er- fiber  oscillator  at  1557  nm 
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was  amplified  in  an  Er/Yb-doped  fiber  amplifier,  then  Raman-shifted  to  1980  nm,  and  finally 
amplified  a  large-mode-area  Tm-doped  fiber  amplifier  [35], 

Two  OP-GaAs  samples  (A3  and  B3)  were  used,  both  0.4-mm-thick  and  3-mm-long.  They  had 
lithographically-defined  QPM  periods  of  1277  pm  (sample  A3)  and  759  pm  (sample  B3).  The 
pump  beam  was  propagating  along  the  [110]  direction  of  GaAs  and  was  polarized  along  [111]  to 
maximize  the  effective  nonlinear  optical  coefficient.  The  focusing  of  the  pump  beam  was 
optimized  to  produce  maximum  THz  power  with  the  resulting  spot  size  (1/e2  intensity  radius)  of 
wo  =  65  pm.  The  generated  THz  beam  was  collimated  with  a  Picarin  plastic  lens  that  had  a  focal 
length  of  50  mm.  THz  power  was  measured  with  a  Si  bolometer  operating  at  liquid  helium 
temperature.  A  black  polyethylene  filter  was  used  to  block  optical  radiation. 

Spectral  properties  of  the  generated  THz  radiation  were  measured  with  a  Michelson 
interferometer.  For  the  sample  A3  the  spectrum  is  centered  at  1.78  THz  and  has  a  width  of  0.3 
THz,  while  for  the  sample  B3  the  spectrum  is  centered  at  2.49  THz  and  has  a  width  of  0.25  THz, 
in  agreement  with  predicted  QPM  peak  position  and  width. 

Figure  5  shows  the  measured  THz  power  versus  the  incident  optical  power  at  1980  nm  for  the 
samples  A3  and  B3,  which  show  very  similar  behavior.  At  the  highest  optical  power  of  2.1  W 
available  at  the  samples  ~  3.3  pW  of  THz  average  power  was  obtained.  The  THz  output  power 
was  quadratic  with  respect  to  the  incident  pump  power  and  did  not  show  any  saturation  effects, 
up  to  the  maximum  pump  power  of  2.1  W  at  the  samples  (peak  intensity  inside  the  samples  1.85 
GW/cm2). 

6.2.  OR  with  fs  pulses  from  an  OP  A/  DFG  system 

Optical  rectification  experiments  were  performed  in  DB-  and  OP-GaAs  using  another 
femtosecond  source,  with  lower  repetition  rate  but  higher  peak  power.  Tunable  pump  pulses 
centered  at  2-4.4  pm,  were  produced  using  a  parametric  amplifier  system  pumped  by  Ti:Sapphire 
pulses  after  a  regenerative  amplifier.  To  achieve  X>  3  pm  wavelengths,  an  additional  difference 
frequency  generation  stage  was  used.  Typical  pulse  durations  were  -100  fs,  repetition  rate  1  kHz, 
and  pulse  energy  up  to  3  pJ  [34],  GaAs  samples  were  not  AR-coated  and  their  main  parameters 
are  listed  in  Table  2. 

The  optical  pump  beam  with  a  beam  size  ranging  between  vv=300  pm  and  1.5  mm, 
propagated  along  the  [110]  direction  (or  its  equivalent)  of  GaAs.  A  Picarin  lens  was  used  to 
collect  the  THz  radiation  to  the  Si  bolometer  (4  K).  To  measure  the  spectral  properties  of  THz 
radiation,  a  Michelson  interferometer  was  used,  composed  of  two  2-inch-diameter  flat  gold 
mirrors  and  a  25-pm-thick  mylar  beamsplitter  (Fig.  6). 


The  power  spectra  of  THz  pulses  were  extracted  in  our  experiment  by  computing  the 
amplitudes  of  the  Fourier  transforms  of  Michelson  interferograms.  Fig.  6  shows  an  original 
interferogram  and  computed  spectra  for  different  samples  and  different  pump  wavelengths.  THz 
spectra  were  noticeably  distorted  by  water  vapor  absorption  lines  (also  shown  in  Fig.6). 
Experimentally  observed  central  frequencies  and  bandwidths  of  THz  pulses  are  in  good 
agreement  with  theoretical  predictions  based  on  known  GaAs  dispersion  relations.  By  changing 
the  pump  wavelength  or  GaAs  QPM  period,  THz  wave  packets  with  central  frequencies  between 
0.9  and  3  THz  were  generated. 
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The  THz  beam  propagated  collinearly  with  respect  to  the  optical  pump  and  was  close  to 
diffraction-limited  (a  pinhole  method  was  used  to  measure  the  far-field  THz  beam  size).  Fig.7 
shows  optical-to-terahertz  conversion  efficiency,  t|thz  ,  for  samples  DB-75  and  A3,  as  a  function 
of  peak  pump  intensity,  /o.  The  pump  beam  size  varied  in  this  case  between  300  and  810  pm. 
One  can  see  that  the  linear  dependence  of  t)thz,  expected  by  theory,  rolls  off  for  I0  >  2  GW/cm2. 
This  roll-off  behavior  was  also  observed  at  similar  intensities  at  shorter,  k~2  pm  pump.  The 
onset  of  saturation  is  most  likely  due  to  nonlinear  refraction  (m)  in  GaAs  which  induces  self¬ 
phase  modulation  and  self- focusing  of  the  optical  pulses.  Indeed,  if  we  take  ri2~  1.5-1  O'4  cm2/GW 
for  GaAs  (Table  1),  at  7o=2  GW/cm2  and  7=6  mm,  the  on-axis  nonlinear  phase  shift  reaches  ~n. 

In  the  sample  DB-75,  0.66  nJ  per  pulse  at  2.2  THz  was  generated  at  w=300pm,  with  2.3  pJ  of 
pump  pulse  energy,  which  corresponds  to  optical-to-terahertz  conversion  efficiency  of  2.9x  10'4, 
internal  efficiency  8.7x  10'4  (the  samples  were  un-coated),  and  internal  photon  efficiency  3.3%. 
These  numbers  are  in  good  accord  with  the  values  calculated  from  Ref.  38. 


7.  Resonantly  enhanced  THz  generation  in  GaAs  using 
DFG 

With  femtosecond  excitation,  the  limit  for  THz  efficiency  was  set  by  parasitic  high-order 
nonlinear  effects,  most  likely  nonlinear  refraction.  The  optical-to-THz  efficiency  can  be 
improved  by  using  a  DFG  scheme  with  longer  (ps)  pulses,  in  the  optimal  case  with  the  duration 
close  to  the  walk-off  time  between  optical  and  THz  pulses.  This  condition,  for  a  10-mm-long 
GaAs  and  pump  at  ~2  pm,  corresponds  to  the  pulse  duration  of  3-4  ps  [38]. 

A  near-degenerate  PPLN  OPO,  synchronously  pumped  by  a  CW  mode-locked  Nd:YV04 
laser  at  1.06  pm  (HighQ  Laser)  was  built  and  investigated.  The  pump  laser  had  repetition  rate  50 
MHz,  pulse  duration  6.5  ps,  average  power  up  to  10W  in  a  TEMoo  mode.  The  PPLN  crystal  (HC 
Photonics),  served  as  an  OPO  gain  medium.  It  was  10-mm-long,  3-mm-wide  and  1-mm-thick, 
had  a  QPM  period  of  14.1  pm,  and  was  AR-coated  for  both  1.06  and  2.1-2.16  pm.  PPLN  was 
designed  for  the  type-II  QPM  three- wave  interaction  (Y-YZ),  so  that  the  signal  and  the  idler  had 
orthogonal  polarizations.  The  high-finesse  standing-wave  OPO  cavity  (Fig.  8)  was  3-m-long  and 
consisted  of  typically  7  mirrors,  all  of  which  were  highly  transmissive  at  1.06  pm  and  highly 
reflective  (>99.9%)  in  the  vicinity  of  the  2.13  pm  degeneracy  wavelength.  Mirrors  M1-M2,  and 
M7  were  flat,  M3-M4  -  had  ROC  of  200mm,  and  M5-M6  -  ROC  of  500  mm.  A  thin-film 
polarizing  beamsplitter  inside  the  cavity  served  as  an  idler-wave  outcoupler,  thus  the  OPO  in  this 
case  was  singly  resonant  (SRO). 

Insets  to  Fig.  8  show  the  signal- wave  far  field  intensity  distribution,  spectral  output  of  the 
OPO  (signal  and  idler)  at  two  different  PPLN  temperatures,  and  OPO  tuning  curve  with  respect 
to  the  PPLN  temperature.  The  OPO  linewidth,  Av~  3cm'1  (100  GHz)  ,  is  preserved  even  when 
the  OPO  crosses  the  degeneracy  point  and  is  close  to  the  one  given  by  the  time-bandwidth  limit. 
In  general,  any  frequency  spacing,  from  0  to  125  cm'1  (0-3.75  THz)  between  the  OPO  signal  and 
the  idler  waves,  can  be  achieved  by  changing  the  PPLN  temperature.  At  the  measured  pump 
beamsize  w=3 1  pm  and  calculated  signal  eigenmode  beamsize  of  w=5 1  pm,  the  OPO  threshold 
was  ~  1  W.  At  the  full  average  pump  power  of  9  W  available  at  the  PPLN  crystal,  the  idler 
output  reached  3  W  with  the  pump  depletion  of  71%.  By  measuring  the  signal- wave  power 


9 


leaking  through  one  of  the  OPO  mirrors,  the  intracavity  resonating  average  signal  power  was 
measured  to  be  as  high  as  120  W.  Thus,  there  was  a  substantial  resonant  -field  enhancement  of 
the  signal  wave,  due  to  low-loss  cavity  (the  total  round-trop  loss  was  estimated  to  be  3%). 


7. 1.  THz  generation  outside  the  OPO  cavity 

In  this  experiment,  the  OPO  was  tuned  near  degeneracy  and  served  as  a  dual-frequency  source  to 
produce  THz  output  via  DFG.  The  signal  and  idler  beams  were  collinearly  focused  into  the  GaAs 
crystal  to  a  beamsize  of  ~  80  pm  and  THz  output  was  measured  with  a  Si  bolometer.  Fig.  9a 
shows  THz  tuning  curves  -  dependence  of  THz  output  vs.  beatnote  frequency  (tuned  by  changing 
PPLN  temperature),  for  4  OP-GaAs  samples  (A  10,  D3,  BIO,  C5  -see  Table  2)  with  different 
lengths  and  QPM  periods.  THz  output  is  maximized  for  each  sample  when  the  THz  frequency 
Q=oo3-(02,  matches  the  QPM  condition  (003  and  (02  are  frequencies  of  the  OPO  signal  and  the  idler 
waves  correspondingly).  Using  a  set  of  few  GaAs  samples,  with  different  QPM  periods,  pretty 
broad  range  of  THz  frequencies  (1-3.3  THz)  can  be  covered.  Also,  one  can  see  that  longer 
samples  give  rise  to  narrower  tuning  peaks,  in  agreement  with  (2). 

Fig.  9b  is  a  log- log  plot  of  THz  average  power  versus  average  pump  power  for  a  5 -mm- long 
OP-GaAs  sample  (QPM  period  704  pm,  center  frequency  2.2  THz).  The  powers  at  the  two  pump 
beams  (OPO  signal  and  idler)  were  slightly  different,  so  geometric  mean  of  the  powers  of  the 
two  beams  was  used  as  the  ‘pump  power’  on  the  abscissa.  The  best  fit  is  a  straight  line  with  a 
slope  of  2,  which  indicates  that  the  THz  output  is  proportional  to  the  product  of  the  powers  of  the 
two  pump  beams.  The  conversion  efficiency  was  about  40%  of  theoretical  [38];  the  discrepancy 
can  be  accounted  for  by  the  beam  clipping,  water  vapor  absorption  in  the  air  at  THz  frequencies, 
and  possible  uncertainty  in  THz  detector  calibration.  Overall,  experiments  with  ps  pulses  confirm 
predictions  of  theory  that  THz  conversion  is  fluence-dependent  and  one  gets  the  same  THz 
efficiency  per  pj  of  pump  pulse  energy  for  both  ps  and  fs  pulses. 

THz  tuning  curves  (THz  peak  frequency  vs.  GaAs  QPM  period)  obtained  in  this  experiment 
are  in  agreement  with  both  OR  results  and  predictions  based  on  known  dispersion  relations  for 
GaAs  (both  optical  and  THz).  Fig.  10a, b  shows  experimental  as  well  as  calculated  tuning  curves 
for  fs  OR  at  4400  nm  pump  and  extra-  and  intracavity  ps  DFG  at  ~  2128  nm  pump. 


7.2.  Intracavity  THz  generation  with  single  resonance 

By  placing  a  GaAs  crystal  inside  the  OPO  cavity,  one  can  take  the  advantage  of  resonant 
enhancement  of  pump  optical  power,  and  thus  improve  optical-to-THz  conversion  efficiency. 
First,  an  experiment  with  a  singly  resonant  OPO  (SRO)  was  performed  with  resonating  ‘signal’ 
wave.  GaAs  crystal  was  placed  near  the  beam  waist  (beamsize  w~200  pm),  formed  by  two 
concave  mirrors  with  500  mm  radius  of  curvature  inside  the  OPO  cavity  (Fig.  11a).  The  signal 

and  idler  beams  propagated  along  the  [1 1 0]  crystalline  direction  of  GaAs  with  polarizations 
aligned  along  [001]  and  [110]  correspondingly.  A  gold-coated  90°-off-axis  parabolic  mirror  (1" 
diameter  and  2"  focal  length),  with  a  3-mm-diameter  hole  in  its  center  was  placed  after  the  GaAs 
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crystal.  It  fully  transmitted  the  optical  waves  and  allowed  -90%  outcoupling  of  the  THz  waves 
due  to  the  large  diffraction  cone  of  the  THz  radiation. 

With  the  AR-coated  6-mm-long  DB-GaAs  crystal  (sample  DB-77,  QPM  period  504  pm),  50 
pW  of  THz  output  was  generated  at  2.8  THz.  THz  average  power  in  this  case  was  measured  with 
a  calibrated  room  temperature  DTGS  detector  from  Bruker  (model  D  201).  Black  polyethylene 
filters  were  used  to  block  the  2-pm  optical  radiation. 

When  the  OC-GaAs  stack  with  comparatively  large  QPM  period,  A,  was  used,  consisting  of 
four  1-mm-thick  (110)  GaAs  plates,  three  peaks  of  THz  efficiency  were  observed.  The  center 
frequencies  of  the  peaks  were  correspondingly  at  0.66,  2.29,  and  3.38  THz  and  their  FWHM 
widths  were  0.39,  0.32,  and  0.18  THz  (Fig.  12).  These  peaks  corresponded  to  the  1-st,  3-rd  and 
5-th  -order  QPM  in  GaAs  with  center  frequencies,  Qn=271c/(  A  An)  xN  (N=  1,3,5  is  the  QPM 
order),  related  approximately  as  1:3:5  if  one  neglects  THz  dispersion.  THz  output  power  as  a 
function  of  THz  frequency  is  shown  on  Fig.  12.  Interestingly,  the  strengths  of  all  three  peaks  are 
comparable:  efficiency  reduction  of  1  IN  due  to  high-order  QPM  process  [52]  is  offset  by  the  Q 
factor  which  scales  as  N2.  Still,  some  difference  in  the  peaks’  strength  can  be  explained  solely  by 
higher  THz  absorption  at  higher  frequencies.  Potentially,  such  3-color  THz  source  can  be  a 
useful  tool  for  spectrally-selective  THz  imaging. 

Figure  13a  shows  the  THz  beam  intensity  profile  after  parabolic  mirror,  reconstructed  from 
scanning  knife-edge  measurements.  Horizontal  and  vertical  1/e2  spot  sizes  of  7.8  and  13.3  mm, 
respectively,  were  close  to  the  diffraction  limit.  THz  beam,  measured  in  the  focus  of  a  f  =50  mm 
Picarin  lens,  was  captured  (Fig.  13b)  by  a  room-temperature  pyroelectric  camera  (Pyrocam  III, 
Spiricon).  The  focused  beam  spot  size  is  2-3  pixels  (1  pixel=100xl00  pm),  which  again 
confirmed  diffraction-limited  performance, 
x 


7.3  Intracavity  THz  generation  with  double  resonance 

In  a  doubly  resonant  OPO  (DRO),  both  optical  waves  are  enhanced  and  the  total  THz  efficiency 
can  be  further  increased.  Schematic  of  the  setup  for  THz  generation  using  QPM  GaAs  inside  a 
DRO  cavity  is  shown  on  Fig.  1  lb.  The  design  of  the  cavity  is  similar  to  that  of  SRO,  except  for 
two,  instead  of  one,  thin-film  polarizing  beamsplitters  and  two  separate  sets  of  end  mirrors  were 
used  to  resonate  simultaneously  the  signal  and  the  idler  waves  [53]. 

To  achieve  double  resonance,  precise  control  of  the  two  cavity  lengths  was  needed.  To 
accomplish  that  goal,  the  ‘dither- and- lock’  technique  was  used,  involving  a  piezo  actuator  at  one 
of  the  mirrors,  a  2-pm  photodetector  and  feedback  electronics.  Interestingly,  when  GaAs  was 
placed  inside  the  DRO  cavity,  self-stabilizing  mechanism  was  observed  [54],  due  to  thermo-optic 
negative  feedback,  similar  to  the  self-stabilizing  observed  earlier  in  a  CW  AgGaS2  DRO  [55].  In 
both  cases  of  active  and  passive  stabilization,  the  DRO  was  able  to  operate  for  >  30  minutes  with 
an  excellent  stability. 

THz- wave  generation  inside  DRO  was  achieved  with  all  three  types  of  QPM  GaAs  samples: 
OP-GaAs  (Z=l()mm,  A=1300pm,  central  frequency  1.27  THz),  OC-GaAs  (L= 8  mm,  A=1060pm, 
1.57  THz),  and  DB-GaAs  (L=6  mm,  A=504pm,  2.8  THz).  In  the  latter  case,  the  average  THz 
power  was  measured  to  be  1  mW,  in  a  diffraction-limited  beam,  which  corresponds  to  better  than 
10'4  optical-to-THz  conversion  efficiency.  The  resonant  intracavity  signal  and  idler  average 
powers  were  measured  to  be  10.2  W  and  17  W  respectively. 
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8.  Demonstration  of  cascaded  optical-to-THz  frequency 
conversion 

The  possibility  of  overcoming  quantum-defect-related  (Manley-Rowe)  limitations  on  the 
efficiency  of  THz  -wave  difference  frequency  generation,  through  cascaded  nonlinear  optical 
interactions  for  single-pass  mixing  of  two  discrete  near-infrared  frequencies,  was  suggested  by 
Cronin-Golomb  [56],  As  the  pump  beams  run  down  the  Stokes  ladder,  the  number  of  terahertz 
photons  should  continually  increase.  A  potential  improvement  by  a  factor  of  5  was  predicted  for 
ZnTe  pumped  at  824  nm  [56],  Cascaded  down-conversion  was  also  predicted  for  the  case  of  OR, 
where  cascading  results  in  continuous  red  shift  of  the  ultrashort  pump  pulse  [38].  This  prediction 
was  experimentally  verified  quite  recently  in  tilted-wavefront  OR  experiments  in  LiNbOs  with 
high-power  fs  pump  pulses  [57,  29]. 

The  spectrum  of  the  optical  radiation  in  the  intracavity  SRO  DFG  experiment  with  the  6-mm 
DB-GaAs  sample  and  center  frequency  2.8  THz  is  shown  in  Fig.  14.  Optical  radiation  was 
extracted  from  the  SRO  cavity  through  one  of  the  highly-reflecting  mirrors,  and  analyzed  using  a 
grating  monochromator.  This  is  the  first  experimental  demonstration  of  the  discrete  optical 
frequency  cascading,  in  the  DFG  experiment,  predicted  in  ref.  56:  in  addition  to  the  signal  and 
the  idler  peaks,  one  can  see  two  additional  red-shifted  peaks,  corresponding  to  the  1-st  and  the  2- 
nd  cascaded  satellites,  Si  and  S2  [54,58],  All  four  peaks  are  spaced  by  the  terahertz  frequency 
(2.8  THz).  Measured  polarization  states  of  the  peaks  are  shown  on  the  top  of  Fig.  14  and  agree 
with  selection  rules,  stating  that  polarizations  of  the  neighboring  peaks  should  be  orthogonal. 

The  number  of  cascading  cycles  depends  on  the  DFG  pump  acceptance  bandwidth  and  on 
THz  frequency.  It  can  be  on  the  order  of  10  for  2-pm  pump,  10-mm-long  GaAs,  and  terahertz 
frequency  of  1-2  THz.  Potentially,  THz  conversion  efficiency  can  significantly  (by  about  a 
factor  of  10)  surpass  the  Manley-Rowe  limit. 


9.  Outlook  for  future  experiments 

The  aspiration  for  our  future  work  is  to  demonstrate  more  than  100%  THz  quantum  conversion 
efficiency  from  IR  to  THz.  To  accomplish  this,  one  needs  to  resonate  not  only  the  OPO  signal 
and  the  idler,  but  the  1-st  and  the  higher  order  Stokes  satellites.  Because  of  the  group  velocity 
dispersion  (GVD),  mainly  due  to  GaAs,  the  OPO  resonator  modes  are  not  equidistant  (Fig.  15a). 
This  leads  to  the  mode  spacing  being  different,  at  the  wings  of  the  wide  (400-500  cm’1)  spectrum 
consisting  of  frequency  clusters  (signal  +  idler  +  Stokes  satellites),  from  that  dictated  by  the 
repetition  rate  of  the  pump  laser,  and  results  in  accumulation  of  extra  phase  making  it  impossible 
to  achieve  resonant  condition.  However,  GVD  can  be  compensated  by  different  methods:  with 
prisms,  by  adding  dielectric  elements  (e.g.  fused  silica)  with  the  opposite  sign  of  the  2-nd  order 
GVD,  and  with  special  GVD  compensating  mirrors.  In  the  case  when  GVD  is  compensated  (Fig. 
15b),  simultaneous  resonance  can  be  achieved  for  the  1-st  and  higher  Stokes  frequencies. 

Fig.  16  shows  longitudinal  frequency  modes  for  the  OPO  signal,  idler,  Stokes  frequencies,  as 
well  as  for  the  THz  radiation.  When  the  modes  of  the  signal,  idler  and  several  (~8)  Stokes 
frequencies  are  resonant,  they  are  expected  to  be  phase-locked,  through  the  action  of  THz  field, 
which  is  common  to  all  optical  frequencies  and  acts  a  phase-locker.  Interestingly,  when  10  or  so 
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optical  frequency  clusters  are  mutually  coherent,  the  optical  pulse,  in  the  time  domain,  is 
expected  to  defragment  into  extremely  short  (~35  fs)  bursts  of  radiation,  within  the  ~5  ps 
envelope  of  the  original  pulse  duration  (Fig.  16b).  Because  of  strong  dispersion  in  GaAs,  this 
mutual  coherence  is  kept  inside  GaAs  for  only  about  a  millimeter  of  length. 


10.  Conclusion 

Periodically-inverted  GaAs  shows  great  promise  for  efficient  generation  of  frequency-tunable 
narrow-bandwidth  (~100  GHz)  THz  wave  packets  via  parametric  down-conversion  in  the  whole 
range  from  0.3  to  4  THz.  Both  femtosecond  (OR)  and  picosecond  (DFG)  optical  pulses  can  be 
used  as  a  pump.  Intracavity  THz  wave  generation  with  resonant  enhancement  of  both  optical 
waves  seems  to  be  the  most  efficient  way  of  producing  THz  output.  When  GaAs  was  placed 
inside  the  cavity  of  a  picosecond  near-degenerate  doubly-resonant  OPO,  lmW  of  average  THz 
output  power  was  achieved.  With  reduced  intracavity  losses  which  can  be  attained  by  lowering 
optical  losses  in  QPM  GaAs,  as  well  as  by  switching  from  the  linear  to  the  ring  cavity  (Fig.  17), 
much  higher  THz  output  can  be  achieved.  Furthermore,  by  allowing  Stokes-shifted  cascaded 
waves  to  resonate  along  with  the  signal  and  the  idler  wave  (which  can  be  achieved  by  special 
design  of  cavity  mirrors  and  by  compensation  of  2-nd  and  3-rd  order  group  velocity  dispersion), 
the  whole  system  can  be  scalable  to  100  mW.  It  was  difficult  to  make  side-by-side  comparison  of 
three  types  of  QPM  GaAs  samples,  because  of  a  very  limited  variety  of  samples.  However,  there 
are  indications  that  they  demonstrate  similar  performance  at  similar  conditions.  Although  OP- 
GaAs  samples  have  insufficient  thickness,  resulting  in  some  THz  beam  clipping,  most  likely  that 
the  OP-GaAs  will  be  work-horse  of  future  applications  because  of  reproducible  technology  and 
precision  of  definition  of  QPM  periods;  for  example,  exotic  gratings  with  aperiodic  of  ‘fanned- 
out’  structures  can  be  easily  created.  OC-GaAs  technique,  on  the  other  hand,  does  not  allow 
creating  a  variety  of  grating  periods,  but  these  samples  are  easy  to  fabricate  and  they  do  not  have 
aperture  limitations. 
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Table  1  Linear,  nonlinear  optical  properties,  and  figures  of  merit  (normalized  to  LiNbO:,)  of 
crystals,  transparent  in  the  0-4  THz  range  and  most  widely  used  for  optical  THz  generation. 


Crystal 

GaAs 

GaP 

ZGP 

GaSe 

ZnTe 

LiNb03 

FiTa03 

Optical  wavelength 
(pm) 

2.1 

1.06 

2.1 

1.06 

0.8 

1.06 

1.06 

Optical  ref.  index 

3.33 

3.11 

3.15 

2.8 

2.85 

2.16 

2.14 

THz  ref.  index 

3.6 

3.31 

3.37 

3.26 

3.2 

5.2 

6.5 

A VI  ^THz~^g  opt 

0.18 

~0 

0.17 

0.34 

~0 

3.0 

4.32 

THz  absorption  at  1-2 
THz,  cm"1 

1 

3.3 

1 

2.5 

9.9 

21.7 

95 

Electro-optic 
coefficient,  pm/V  [43] 

r41=1.5 

r4i=0.94 

r41=1.6 

r22-1.58 

r4i=4 

r33=28 

r33=27.7 

Nonlinear  coefficient, 
pm/V 

d14=46.1 

d14=21.7 

d36=39.4 

d22=24.3[14] 

p* 

II 

On 

o\ 

d33=152.4 

d33=145.2 

Nonlinear  refractive 

150  [44] 

20  [45] 

40* 

45* 

120[46]@1.06 

0.91  [48] 

0.37* 

index  n2,  10"15  cm2/W 
@  X  (pm) 

@2.1 

@0.78 

@2.1 

@  1.06 

71[47]@0.8 

@1.06 

@1.06 

FOMi,  long  pulses 

0.83 

0.06 

0.68 

0.13 

0.25 

1 

0.21 

FOM2,  ultrashort  pulses 

0.64 

1.67 

0.55 

0.13 

6.04 

1 

0.64 

FOM3,  Kerr-limited 
conversion 

0.014 

0.005 

0.047 

0.004 

0.0037 

1 

0.416 

*  obtained  indirectly  through  the  bandgap  scaling  law  of  ref.  48 


Table  2  Parameters  of  the  QPM  GaAs  samples  used  in  experiments. 


Sample 

AR  coated 

QPM  type 

Aperture  (mm2) 

Fength 

(mm) 

QPM  period 
(pm) 

DB-75 

no 

DB-GaAs 

10x10 

6 

504 

DB-77 

yes 

DB-GaAs 

10x10 

6 

504 

A3 

no 

OP-GaAs 

0.4x3 

3 

1277 

A10 

no 

OP  -GaAs 

0.4x3 

10 

1277 

B3 

no 

OP  -GaAs 

0.4x3 

3 

759 

B5 

no 

OP  -GaAs 

0.4x3 

5 

759 

B10 

no 

OP  -GaAs 

0.4x3 

10 

759 

C5 

no 

OP  -GaAs 

0.4x3 

5 

564 

D3 

no 

OP  -GaAs 

0.4x3 

3 

932 

3-07-1300- 

10A 

yes 

OP  -GaAs 

0.8x5 

10 

1300 

OC-4-lOOO 

yes 

OC  -GaAs 

50x50 

4 

2000 

OC-15-530 

yes 

OC  -GaAs 

5x5 

8 

1060 
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Figure  1  Schematic  illustration  of  collinear  THz-wave  generation  in  a  nonlinear  crystal  with 
periodically  inverted  sign  of  x(2)-  (a)  Optical  rectification  with  femtosecond  pulses,  (b)  difference 
frequency  generation  with  two  picosecond  pulses  (Q=CD3-oe>2). 
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Figure  2  THz  power  absorption  coefficients  for  GaAs  [2],  GaP  [15],  ZGP  [39],  GaSe  [14],  ZnTe  [40], 
LiNb03  [41],  and  I  .i'l  aO;  [2],  measured  at  room  temperature. 
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Figure  3  Nonlinear  coefficient  du  for  GaAs,  as  a  function  of  THz  frequency  a>i,  provided  that  the  two 
pump  frequencies  (co2,  cd3)  are  larger  than  the  optical  phonon  frequency  [17]. 
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Figure  4  OP-GaAs  samples  for  the  optical  THz  generation  with  Z=10mm  and  the  QPM  period  of  1300 
pm.  (a)  Top  view,  (b)  Transmission  picture  of  the  cross  section  obtained  with  a  near-infrared  (k~  1  pm) 
microscope,  from  the  polished  face,  as  viewed  by  the  propagating  optical  beam;  (c)  transmission  picture 
of  the  grounded  side  face;  one  can  see  boundaries  between  inverted  domains.  The  usable  thickness  of  the 
samples  in  this  case  is  ~0.8  mm. 
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Figure  5  Schematic  of  the  OR  THz  experiment  with  a  2 -jam  fiber  laser.  THz  average  power  vs.  laser 
pump  power  for  the  OP-GaAs  samples  A3  (filled  circles)  and  B3  (open  circles).  Straight  line  is  the  best  fit 
for  the  A3  data. 
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Si  bolometer  (4K) 


Figure  6  (a)  Experimental  setup  for  Michelson  interferometry  and  (b)  typical  Michelson  interferogram. 
Power  spectra  of  THz  pulses  extracted  from  interferograms,  for  different  samples  and  pump  wavelengths: 
(c)  BIO,  pump  at  4.4  pm,  (d)  DB-75,  pump  at  4.4  pm,  and  (e)  DB-75,  pump  at  2  pm.  Shown  in  gray  is 
transmission  of  a  20  cm  path  of  standard  air  taken  from  the  HITRAN  database. 
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Figure  7  Optical-to-terahertz  conversion  efficiency  as  a  function  of  peak  pump  intensity  for  the  sample 
DB-75  (2.2  THz,  circles)  and  A3  (1.5  THz,  triangles).  The  pump  beam  size  was  810  pm  (open  circles), 
520  pm  (closed  circles),  590  pm  (triangles),  and  300  pm  (crossed  circles).  The  pump  wavelength  was  3.5 
pm  (open  circles  and  triangles)  and  4.4  pm  (filled  and  crossed  circles).  Dashed  lines  -  linear  fits. 
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Figure  8  Near  -degenerate  type-II  PPLN  OPO,  synchronously  pumped  by  ps  pulses  at  1.06  pm.  Insets: 
‘signal’  wave  far  field  intensity  distribution,  spectral  output  of  the  OPO  at  two  different  PPLN 
temperatures,  and  OPO  tuning  curve  with  respect  to  PPLN  temperature. 
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Figure  9  (a)  THz  output  vs.  frequency  in  the  DFG  experiment  outside  the  OPO  cavity,  for  4  OP-GaAs 
samples  with  different  lengths  (3-10  mm)  and  QPM  periods  (564-1277  pm),  (b)  THz  average  output 
power  versus  incoming  average  pump  power  for  a  5-mm-long  uncoated  OP-GaAs  sample  (QPM  period 
704  pm,  center  frequency  2.2  THz).  Straight  line  is  the  best  fit. 
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Figure  10  THz  peak  frequency  vs.  GaAs  QPM  period  for  (a)  fs  OR  at  4400  nm  pump  and  (b)  ps  DFG  at 
2128  nm  pump.  Dashed  lines  -  calculated  from  known  GaAs  dispersion. 
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Figure  11  Setup  for  THz  generation  using  GaAs  inside  the  cavity  of  a  near  -degenerate  synchronously- 
pumped  PPLN  OPO.  (a)  Singly-resonant  OPO  (SRO),  (b)  doubly-resonant  OPO  (DRO).  Off-axis 
parabolic  mirror  with  a  hole  in  the  center  was  used  to  extract  the  THz  output. 
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Figure  12  THz  output  as  a  function  of  frequency,  in  the  case  of  intracavity  THz  generation  using  an  OC- 
GaAs  stack  consisting  of  four  1-mm  (1 10)  optically-contacted  plates.  The  peaks  correspond  to  the  1-st,  3- 
rd  and  5-th  -order  QPM  in  GaAs. 
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Figure  13  (a)  Collimated  THz-beam  intensity  profile  reconstructed  from  scanning  knife-edge 
measurements,  (b)  Focused  (f=50  mm  lens)  THz-beam  intensity  profile  measured  by  a  room-temperature 
pyroelectric  camera  (1  pixel=100  xlOO  pm). 
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Figure  14  Spectrum  of  the  optical  radiation  (not  to  scale)  in  the  intracavity  SRO  DFG  experiment  with 
the  6-mm  DB-GaAs  sample  (center  frequency  2.8  THz).  In  addition  to  the  signal  and  the  idler  peaks,  one 
can  see  two  additional  red-shifted  peaks,  corresponding  to  the  1-st  and  the  2-nd  cascaded  satellites  Si  and 
S2.  Polarizations  states  of  the  peaks  are  shown  on  the  top. 
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Figure  15  (a)  Case  of  non-zero  group  velocity  dispersion  (GVD).  The  OPO  resonator  modes  are  not 
equidistant  and  it  is  impossible  to  achieve  resonant  condition  for  the  1-st  and  higher  Stokes  frequencies, 
(b)  When  GVD  is  eliminated,  simultaneous  resonance  can  be  achieved  in  the  whole  frequency  range. 
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Figure  16.  (a)  Frequency  modes  for  the  OPO  signal,  idler,  Stokes  frequencies,  as  well  as  for  THz  wave, 
(b)  Time  domain  picture.  When  the  modes  of  the  signal,  idler  and  several  (-8)  Stokes  frequencies  are 
resonant  and  phase-locked,  the  optical  pulse  consists  of  extremely  short  (-35  fs)  bursts  of  radiation. 
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Figure  17.  Conceptual  schematic  of  the  setup  for  resonantly-enhanced  THz-wave  generation  in  a  ring 
OPO  cavity. 
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